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ON THE STABILITY OF THREE-DIMENSIONAL BOUNDARY LAYERS 
PART 1: LINEAR AND NONLINEAR STABILITY * 

ERIK JANKE AND PONNAMPALAM BALAKUMAR 


Abstract. The primary stability of incompressible three-dimensional boundary layers is investigated 
using the Parabolized Stability Equations (PSE). We compute the evolution of stationary and travelmg 
disturbances in the linear and nonlinear region prior to transition. As model problems, we choose Swept 
Hiemenz Flow and the DLR Transition Experiment. The primary stability results for Swept Hiemenz Flow 
agree very well with computations by Malik et al. For the DLR Experiment, the mean flow profiles are 
obtained by solving the boundary layer equations for the measured pressure distribution. Both linear and 
nonlinear results show very good agreement with the experiment. 
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Subject classification. Fluid Dynamics 

1. Introduction. The transition process without bypass in a three-dimensional boundary layer flow 
can be described by five main phases. First, in the receptivity phase, disturbances determined by the outer 
conditions (surface roughness, freestream turbulence, acoustic noise) are entering the boundary layer. In the 
region of their introduction, a wide spectrum of disturbances is present. Many of these initial disturbances 
decay, however, and only a few are amplified in the downstream flow. Second, the phase of exponential 
growth sees a slow amplitude growth of the few unstable modes. Due to its linearity, this phase can be well 
described by following the most unstable mode. Third, in the phase of nonlinear interaction, the disturbance 
amplitudes are large enough to interact nonlinearly, the uniform spanwise mean flow is modulated by the 
disturbances, and a disturbance saturation into an equilibrium stage can be observed. The location of the 
saturation onset, as well as the saturation amplitude level are strong functions of the initial conditions of 
disturbance amplitude, frequency and phase, which are determined by the receptivity mechanism. In this 
third phase, the mean flow profiles are being distorted, the well known co-rotating crossflow vortices are 
fully developed, and the boundary layer thickness varies strongly in spanwise direction. Due to the strong 
distortion of the boundary layer, inflectional mean flow profiles develop and we can distinguish a fourth 
phase, in which the boundary layer becomes unstable to three-dimensional high-frequency disturbances. 
The frequencies observed in this phase are an order-of-magnitude higher than the primarily unstable ones 
and this mechanism is usually referred to as secondary instability. Finally, an explosive growth of these 
high-frequency modes initiates the fifth phase, the breakdown to turbulence. 

In the research of purely crossflow instability dominated incompressible boundary layers, significant 
progress was made within the last years. Pioneering experiments, in particular performed in Gottingen, 
Germany and in Tempe, Arizona, U.S.A., are accompanied by computational work that fives more and more 
up to the ch alle nge of guiding the experiments. Using the Parabolized Stability Equations (PSE), Malik 
et al. [17] showed for Swept Hiemenz Flow that including traveling disturbances with a small amplitude 
in the analysis can significantly decrease the saturation amplitude and influence the onset of saturation in 
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comparison with a purely stationary computation. Incorporating the receptivity to surface roughness in 
his PSE computations for the DLR Transition Experiment and also introducing initially small amplitude 
traveling modes, Bertolotti [4] arrived at results that are in excellent agreement with the measurements. 
Joslin and Streett [13] showed for stationary crossflow vortices that prior to the nonlinear interaction of 
the disturbances, linear superposition of the smaller amplitude modes can explain wave angle and phase 
adjustments observed in experiments by Arnal and Juillen [1].' Also, results by Kohama et al. [14] and 
Reibert et al. [18] for a swept wing boundary layer in very benign environmental conditions showed the 
existence of a high-frequency secondary instability that is caused by the inflectional deformation of the mean 
flow profiles. 

In the present paper, we study the linear, nonlinear and secondary stability of the Swept Hiemenz Flow 
model problem and of the crossflow dominated boundary layer flow in the DLR Transition Experiment. For 
the sake of a more systematic presentation of the results, we divide the paper in two parts. Starting from 
a mean flow computation and a nonlinear stability analysis, we obtain the equilibrium solutions for each 
problem using the PSE in Part 1. In Section two of Part 1, we describe the present formulation for the mean 
flow computation and the primary stability analysis using the PSE. Sections three and four contain results 
from linear and nonlinear PSE computations for Swept Hiemenz Flow and the DLR Transition Experiment, 
respectively. We then determine the local stability characteristics in the equilibrium region using temporal 
Floquet analysis in Part 2 of this paper. There, the main focus is on investigating the existence of multiple 
roots in the eigenvalue spectrum of the modified mean flow in the region of nonlinear disturbance saturation. 

2. Formulation. In this section, we present the mean flow computation and describe the PSE formula- 
tion. Figure 2.1 shows the coordinate systems that we use in the computations. The coordinates (xj ,x 2 ,x s ) 
represent the body-fixed coordinate system, and the (xj , x 2 , x$) coordinates stand for the vortex-oriented 
system which is applied in the secondary stability analysis. 



Fig. 2.1. Body-fixed (xj } X2,X3 ) and Galilean coordinate system (xj , £$ ) 


2.1. Basic Flow. For a general quasi-three dimensional geometry and a corresponding coordinate sys- 
tem (Figure 2.1), we define a steady incompressible basic flow Q 0 (x, ,x s ,z s ) = { U 0 , W 0 , V 0 ,P 0 } T - Assum- 
ing no variations of the basic flow in the spanwise direction, introducing a stream function § (Equation 2.1), 
and defining the mean flow components according to Equation 2.2: 


( 2 . 1 ) 

( 2 . 2 ) 


= v/W; U;f{x’,,r)) ; T) = 


; m 


dr: 


Vtrxfjuj ’ u* dx; 


U'o = dtjdx ■; ; v; = -d*jdx\ ; w* 0 = ■ g {x-, v ) , 
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we can write the boundary layer equations and the boundary conditions for a quasi-three-dimensional flow 
in the following form, where primes denote differentiation with respect to rj : 


r + ■ r + ■» • u -/") = *>•(/'• - f" 


df 


dx t ’ ’ 


(2.3) 


„ m + 1 , dg , df 

9 +—I-9 =*<•(/ ’ - 9 - aS7> ’ 


f =f = g = 0 at 1 1=0 ; /' = g = 1 at r/->oo 


Here, U* is the inviscid velocity in the chordwise direction, is the constant inviscid velocity in the 
spanwise direction, and m is the similarity parameter. This set of equations is solved using the two-point, 
fourth-order-accurate compact scheme ([16]) in the wall normal direction and a second-order-accurate three 
point upwind scheme in the streamwise direction. 

2.2. Parabolized Stability Equations. In the Parabolized Stability Equations approach, one at- 
tempts to construct an approximate solution of the Navier-Stokes Equations. The concept was introduced 
by Herbert and Bertolotti [10] and applied to linear and nonlinear Blasius boundary layer flow by Bertolotti 
et al. [3]. Now, it has been well developed and applied to two- and three-dimensional incompressible and 
compressible boundary layer flows [6, 2, 17]. 

If we write the total flow quantity Q={V, P} T ={ U, W, V,P} T as the summation of the mean flow 
Qo={V 0 , Po } T ={ U 0 , W 0 , Vo , Po } T and the disturbance q t , Pi } T ={ui , Wj , v t , pi } T , the Navier-Stokes 
equations for the disturbances quantities become : 


(2.4) 


V * = 0 , 


(2.5) + (vj ■ V)v t + (Vo ■ V)vj + (vj ■ V)V 0 = - Vpi + . 

The central idea of the PSE formulation is to split the disturbance quantities into a slowly varying shape 
function and an oscillatory wave part as: 


(2.6) 


Ql ( j 1 0 — 


E 


Q.1 ,mn 


( \ * f Of ,t»n dxi +im(3i xt —tnurt 

(xi,xs)e J 


n=—oo m= — oo 


where j3i is the spanwise wave number, u; is the frequency, ai is the wave number in the streamwise 
direction, and qt , m n is the vector of the shape functions. The indices m and n are the Fourier summation 
indices and also denote the mode number (m,Ti). Substituting this into Equations 2.4, 2.5 and neglecting 
the second derivative of the shape function in the Xt -direction, we obtain the PSE (Equation 2.7) for each 
Fourier component. In Equation 2.7, the vector S mn represents the Fourier component of the nonlinear term 
(vx • V)«i, and A, B, C, D, E are (4 x 4) coefficient matrices. 


(2.7) 


d 2 qi,r 


dx§ 


+ B, , 


dq i,r, 


4“ C mn qt , rnri — B r 


dqt )f i 


+ B-n 


d 2 qi,v 


+ Sri 


dx 3 ' - - dxt ' T, ‘" dudxs 

We note that the linear PSE are obtained by setting the indices m and n in Equation 2.7 to one, and by 
dropping the nonlinear term S mn . Further, assuming that the amplitude part q t of the disturbance quantity 
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is no longer a function of the streamwise location, but depends only on the normal coordinate, the two 
remaining terms on the right hand side of Equation 2.7 vanish and we obtain the linear stability equation. 

In the spatial framework of the PSE, we have the unknown wave numbers oil , 0 j , and the frequency 
^ bi addition to the unknown shape functions Qj . Pixing the real frequency and the real spanwise wave 
number 0 , , a relation to determine the complex streamwise wave number a, is still missing. This forms 
the main difficulty in the PSE formulation, since both the shape functions and the phase of the disturbance 
quantities depend on the streamwise coordinate (Equation 2.6). In a non-parallel mean flow, the different 
physical quantities grow at different rates, and thus, one can determine the growth rate crj =-ai,imag only 
from the computed wave number aj based on some quantity (e.g., velocities, pressure, or energy). Usually, 
the wave number oc i is computed at the location in the boundary layer where the disturbance quantities 
become maximal. This location slowly varies in the marching direction. The procedure to compute q, ata 
streamwise station is described as follows. Starting with an initial guess that is found by solving the linear 
eigenvalue problem at x, = x 0 , one marches to the next station x t = x 0 + A xj and solves for the shape 
functions qi assuming that at j (x; = xo + Axj) = aj (xo). Approximating the change in the shape functions 
with a Taylor Series at xj = xq + Ax, one can derive Equation 2.8 to update a , at x t = x 0 4- Ax t as : 


( 2 . 8 ) 


a l ,ne w — ,old + 


(dqA 

■ q t \dxj ) 




where qj stands for any physical quantity in the flow field. Updating the shape functions and iterating until 
the change in a t is less than a defined tolerance, the solution at this station is obtained and the procedure 
repeated at the next streamwise station. The PSE method is very efficient, since it takes only a few iterations 
°n oli to obtain an accurate solution, provided that the gradients in the flowfield are moderate. In our formu- 
lation, the PSE are solved by introducing a new vector of unknowns q t ={iij , dui /dx 3 , u 2 , du 2 / dx 3 ,u 3) p t } T 
that allows for a discretization of the wall-normal direction using the two-point, fourth-order-accurate com- 
pact scheme [16]. For the derivatives in the streamwise direction, we employ a first-order-accurate two-point 
upwind scheme. 


3. Swept Hiemenz Flow. 

3.1. The Basic Flow. The model problem of Swept Hiemenz Flow represents an exact solution to the 
Navier-Stokes Equations and can be used to model the leading edge region of a swept wing. The inviscid 
velocity distributions are given by: 

(3- 1 ) U' = c-x ; ; U* = - C . a; ;; W = = const. . 


Further, we define the nondimensional frequency F, the constant local length scale /*, and the stream- and 
spanwise Reynolds numbers Re and Re, respectively, in Equation 3.2. The mean flow profiles can be obtained 

(3.2) F = 2nvT/(W^) 2 ; ^ ; Re = ; Re = 

l/* v* 

by solving the Falkner-Skan-Cooke equations ([15, 17]), which are the homogeneous version of Equations 2.3 
with m=l. The results are obtained for 7?e=500, F= 0.7510 -4 and a spanwise wave number of 0j —0.4. 
The chosen value for F is close to the most amplified frequency for this flow. 

3.2. Nonlinear Analysis. Following the previous work by Malik et al. [17], we investigate the evolution 
of stationary and traveling disturbances with and without wave interaction. In particular, we consider the 
interaction of a stationary and a traveling disturbance of similar initial amplitude, as well as the interaction of 
a stationary disturbance with larger amplitude and a traveling disturbance with smaller amplitude. Figure 3.1 
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Fig. 3.1. Growth rates based on disturbance- energy 


shows the growth rates for different initial amplitude levels A. The rms- amplitude A is defined by A = 
y/{uj + wf )max/ Woo, and the computed initial disturbances are introduced at Re~ 186. The computations 
are performed including eight modes in the spanwise and frequency domain and the results agree with the 
results by Malik et al. [17]. It can be seen from Figure 3.1 that the interaction of a larger amplitude stationary 
with a smaller traveling disturbance leads to an earlier saturation of the stationary disturbance. On the other 
hand, the nonlinear evolution of the traveling disturbance is not affected by the stationary vortex when both 
initial modes are of similar amplitude. 

The influence of including a small amplitude traveling disturbance in the computations is presented in 
Figures 3.2 and 3.3, where we show the evolution of the stationary and traveling disturbance amplitudes. 
In Figure 3.2, we plot the results for the purely steady case, and in Figure 3.3 for a case where a small 
amplitude traveling disturbance interacts with the stationary vortex. The latter case will be referred to as 
the “lower-frequency” case from now on. For the purely stationary case, we observe that both the - and 
Wi -components are fully saturated at Re ~ 550, and that the saturation amplitude levels are A s — 24% for 
the Ui -component and A s — 17% for the wi -component. For the “lower-frequency case with a primary 
frequency of F-0.75-10 -4 , we see that the stationary disturbances saturate at Re ^ 480 and the amplitude 
levels are A, ~ 10 % and A 8 ~ 8 % for the u t and w t -component, respectively. Also, since the traveling 
disturbances grow longer in that case (see curve 4a in Figure 3.1), they reach larger amplitude levels than 
the stationary disturbances. 

Next, we document various velocity profiles for the purely stationary case. At the streamwise location 
where the stationary disturbances are saturated (f£e=546), we plot the shapes of the individual Fourier 
modes for the velocity components along the stationary vortex (u/i), and perpendicular to it («i) 111 the 
Galilean coordinate system (xj , X 2 , £ 3 ) in Figure 3.4. Obviously, the first seven modes of the component 
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along the vortex contain important information, and thus, it is advisable to include at least eight modes in 
the computation. We also observe that the fundamental mode (3i along the stationary vortex is very large 
at this Reynolds number (A ~ 30%) and shows the typical structure observed in the experiments just before 
the secondary instability mechanism causes transition. 



Fig. 3.2. Amplitudes of the primary disturbances for the purely stationary case 



Fig. 3.3. Primary disturbance amplitudes for the “lower- frequency case ” (solid : traveling mode, hollow : steady mode) 
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Fig. 3.4. Shape functions of individual modes at Re- 546 

Figure 3.5 shows the total velocity component tangential to the stationary vortex W 2 in the xj , x 5 -plane 
at four different locations within one wavelength perpendicular to the vortex. Finally, in Figures 3.6a)- 
d), we show the toted velocity components (mean flow and disturbance) along the stationary vortex {W 2 ) 
and in the wall-normal direction ( V 2 ) at two streamwise positions and at different normal locations. At 
Re— 486, a developed peak- valley structure of the W 2 -component is observed, and spanwise gradients in 
the normal velocity component are seen as well. At Re— 546, however, an intricate spanwise variation of 
both plotted velocity components is noted. The increased inflectional character of the spanwise velocity 
profiles corresponds to the growing instability of the flow. Further, the strong spanwise gradients in the 
wall-normal velocity component (see Figure 3.6(d)) indicate the presence of a strong localized fluid-exchange 
mechanism between the wall region and the outer flow. As in Figure 3.5, the inflectional character of the 
profiles is clearly observed. These strongly inflectional profiles in all coordinate directions are the origin of 
the secondary instability that will be studied in Part 2 of this paper. 



Fig. 3.5. Total flow quantity at different spanwise stations for Re- 546 
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a) 


b) 




Fig. 3.6. Horizontal cuts of the modified mean flow at two itreamwise location » ( a), b) at Re - 4 86; c), d) at Re=5f6 ) 

4. DLR Transition Experiment. In this section, we investigate the linear and nonlinear evolution of 
stationary and traveling disturbances for the experiment reported by Deyhle and Bippes [9]. The experimen- 
tal pressure gradient across the entire swept flat plate is favorable, and hence, the transition in the boundary 
layer is dominated by a purely crossflow type instability. Because of the wide variety of the investigated 
parameters and the well documented experimental data, this transition experiment is of high value for a 
numerical investigation of the crossflow instability mechanisms. 

4.1. The Basic Flow. The experiment was performed for a freestream velocity of <?£,=19 m/s, a 
sweep angle of 0=45°, and a chord length of the model of c*=0.5m. In order to compensate for blockage 
effects and problems with the simulation of infinite swept wing conditions using end plates, Deyhle and 
Bippes [9] suggest the use of (?^ o =20. 5 m/s and 0=43.5° instead. Our computations are done for the latter 
parameters and a Reynolds number based on the chord length of Re^ =683,000, where Re x = Q' x ■ c’ jv* . 
Further, we use U e 0 and = • x' I0 j 11*, 0 ) 1 ' 2 as the velocity scale and length scale, respectively, where 

Xj 0 and U* 0 are the streamwise coordinate and the freestream velocity in the x, -direction at the location 
where we introduced the disturbances. The turbulence level during the experiment is specified as Tu=0.15%. 
Even though this value represents a fairly benign environment, it is still high in comparison to the turbulence 
level present in the experiments performed by Reibert et al. [18] and Saric et al. [19] where they quote a 
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value of Tu=0.04% at a speed of Q^=20 m/s ([7]). In the DLR Experiment, traveling modes are present 
in the flow and one needs to consider these in the computations when comparing with the experimental 
results. Further, the fact of the initial presence of traveling modes complicates the isolation of the physical 
mechanisms considerably. 

In contrast to previous work which used a Falkner-Skan-Cooke (FSC) similarity solution to fit the 
measured data, we compute the mean flow profiles directly from the measured outer velocity distribution 
that was obtained by courtesy of the DLR (Bertolotti [5]) and is given in Figure 4.1. In Figure 4.2, we 
compare the boundary layer profiles along the inviscid streamline Ut and the crossflow velocity profiles U c 
resulting from the present approach and the FSC solution (Deyhle and Bippes [9]) with the experiment. We 
note a good agreement of the computed with the measured profiles. 



FlG. 4.1. Measured outer velocity distribution and angle of inviscid streamline 



Fig. 4.2. Mean flow profiles tangential and perpendicular to the streamline at xj/c-O.J 
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4.2. Linear Stability Analysis. First, we conduct a local stability analysis of the undisturbed mean 
flow. In the spatial formulation of the local stability analysis, we consider a normal mode of the form : 

(4-1) qi(xi ,x s ,x 3 ,t) =q 1 (x 3 )e ,a,x,+l0,x, '~ iwt , 

where a, is the complex wave number in the chordwise direction, 0, is the real wave number in the spanwise 
direction, and u> is the real frequency. The wave vector k^ai, the wavelength A and the wave angle are 
defined in Equation 4.2. 

(4-2) kreal = {« 1 ,rcah 01 } T ! A = — j= ■ = = — ; = tan ' ~ 1 — — 

\/ a i, real +0*1 ' k ™‘' Ql ’ reo ' 

According to Mack (1984), two famihes of most unstable modes are present in three-dimensioneil boundary 
layers. Figure 4.3 depicts the propagation direction of the two unstable famihes. It also defines S' as the 
angle between the wave vector and the inviscid streamline. 



Fig. 4.3. Definition of the wave angle and the angle of the inviscid streamline & 

In order to compare our linear results with experimental data, we consider the experiment performed 
by Deyhle et al. [8]. First, we show two unstable disturbance famihes for the arbitrary frequency of /*=150 
Hz versus the wavelength in the spanwise direction A*,. From Figure 4.4, it is seen that the more unstable 
disturbances (/?, > 0) propagate in the direction of the centrifugal force induced by the inviscid streamline. 



Fig. 4.4. Two unstable disturbance families (zijc=0.1, f* =150 Hz) 
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Second, we compute wavelength, wave angle and phase speed at a streamwise location of £j/c 0.7 as 
functions of the frequency by solving the spatial eigenvalue problem. Considering the temporal eigenvalue 
problem, the group velocity was determined. In Figures 4.5 and 4.6, we plot the wave angle & relative to 
the angle of the inviscid streamline, the wavelength A = 27r/\k rea i\ , the phase speed C p hase — v/\Keai\ 
the group velocity C gro up = \du/dk r€a i\ for both families. The results are presented for the most unstable 
wave at a fixed frequency. Between the two modes with positive and negative (3i , the mode with a positive 
spanwise wave number has the largest growth rate. It is this family with a positive f} t that seems to match 
the experimental results reasonably well. Considering the results for the wave angle, the wave length and 
the phase velocity, we note a better agreement with the experiment than for the group velocity results. 
However, Deyhle et al. [8] report experimented uncertainties in determining the group velocity, which might 
also account for the deviations in these results. 



Fig. 4.5. Experimental results and local theory at zj j c — 0.7 



Fig. 4.6. Experimental results and local theory at zi/c — 0.7 
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Next, and reconsidering the experiment described in Section 4.1, we show the results of a local stability 
analysis performed at different streamwise locations. Thereby, the critical points and the locally most 
amplified steady and traveling disturbances are determined. 

In Figure 4.7, the nondimensional growth rates o j for stationary disturbances are plotted versus the span- 
wise wave number /3j at different chordwise positions. It can be seen that the local theory predicts an increase 
of the unstable spanwise wave number range in the downstream direction for the stationary disturbances. 
The critical point for the stationary disturbances is found at x, /c=0.067 for =0.174 (A* t = 10.5mm). The 
locally most amplified wave number is determined as /3, =0.295 at x, /c=0.40, where the chordwise wave 
number is aj =(-0.33509,-0.01401). Table 1 contains the locally most amplified spanwise wave numbers, as 
well as the corresponding growth rates, dimensional spanwise wavelengths and wave angles. 



Fig. 4.7. Stationary growth rates from a local analysis 


Table 1: Locally most amplified eigenpairs at different chordwise locations (/*=0 Hz) 


Xj/c 

(j i = max . 

,cr / =max [wi?7l] 

ftl y( rj — max. 

dr 

* cri — max. 

0.07 

0.00033 

10.43 

0.177 

87.60° 

0.10 

0.00289 

10.92 

0.193 

87.45° 

0,20 

0.00858 

11.46 

0.234 

86.89° 

0.30 

0.01209 

11.31 

0.269 

86.28° 

0.40 

0.01401 

11.23 

0.295 

86.09° 

0.50 

0.01387 

11.57 

0.308 

85.70° 

0.60 

0.01310 

12.27 

0.309 

86.54° 

0.70 

0.01279 

13.02 

0.308 

86.48° 

0.80 

0.01317 

13.32 

0.313 

86.18° 


It is observed that the wave angle of the most amplified stationary disturbances remains almost constant 
in the chordwise direction. The constant phase lines of the most unstable stationary disturbances are inclined 
at angles less than five degrees with respect to the inviscid streamline. 

Third, we introduce stationary and traveling disturbances separately at x t /c=0.06 and perform a linear 
PSE analysis to determine both the most amplified frequency and spanwise wave number. In Figure 4.8, 
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we show a comparison of the linear PSE computation for the stationary disturbance components with 
measured data of Deyhle and Bippes [9]. The results are presented for the velocity components in the 
directions tangential and perpendicular to the in viscid streamline. Both the locally parallel and linear PSE 
computation capture the essential features of the measured disturbance profiles like the triple-peak structure 
and the location of the crossflow-component maximum very well. Despite the good qualitative agreement of 
the computed crossflow component with the experimental results, we note that the magnitudes obtained from 
measurement and computation differ by a factor of two. This was also reported by Deyhle and Bippes [9] 
and is attributed to the presence of small spanwise gradients in the mean flow. Further, the local theory is 
predicting a slightly too high peak location of the tangential component, whereas the linear PSE matches 
the experimental data better. 



Fig. 4.8. Stationary disturbance profiles from experiment and computation at xj/c-0.6 

Finally, Figure 4.9 shows both the calculated linear and nonlinear growth rates based on the u t - 
component, as well as the ^-factors for the most amplified modes. The obtained values for the most amplified 
spanwise wavelength A* f =12mm, the most amplified frequency /*=178Hz, and the maximum JV-factor of 
JV~16 for the traveling disturbance agree exactly with previous work [4]. We also note from Figure 4.9, that 
the traveling disturbances are much more amplified than the stationary disturbances. 

4.3. Nonlinear Analysis. Since our formulation does not include the receptivity to freestream tur- 
bulence or surface roughness, we match the initial amplitudes of the disturbances with those given by 
Bertolotti [4] at x, /c=0.1. The amplitude is defined as A = | m ai/ ^oo- F° r — 0.04% in the body-fixed 
coordinate system, we obtain the same value for the disturbance component along the inviscid streamline 
Utang,avg = )/* *t x 1 /c=Q.l as shown by Bertolotti. Here, the subscripts max and 

min stand for the maximum and minimum value of an average in the spanwise direction. For the nonlinear 
computations in this section, we include eight modes in spanwise direction and four modes for the frequency. 
As reported by Bertolotti [4] and in Section 3.2, we find that the purely stationary disturbance saturates at 
a higher amplitude level and a later streamwise location than in the experiment. Therefore, we include a 
small amplitude traveling disturbance with the most amplified frequency of /* — 178Hz in the computations 
(A t , u = 0.004%). 
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Fig. 4.9. Linear and nonlinear growth rates and N -Factors for =12mm , /* -178Hz 



In the nonlinear cases, the wave number and frequency domains are modeled using 8x4 modes for /3 2 
and ui, respectively. Allowing for the nonlinear interaction of these disturbances, the saturation onset is 
shifted upstream and the saturation amplitude decreases. Figure 4.10 shows the measured and computed 
evolution of the disturbances along the chord. We observe that the saturation onset is predicted slightly 
late, but the experimental saturation amplitude is matched quite well. 

Next, we describe the spatial development of the crossflow vortices. Figure 4.11 shows the u, -disturbance 
component at xj j c— 0.5, 0.6, and 0.7 for both the purely stationary and the interaction case together with 
the experimental data at x , j c=0.6. It is seen that the purely stationary vortices from the computation at 
xj / c= 0.6 are further developed than the measured vortices. Both their extension in the normal direction 
and their degree of distortion are overpredicted. In the interaction case, however, the normal extension of 
the vortices has decreased to the experimental value, but the structure of the vortices is more distorted 
than shown in the experimental results. This difference may be attributed to the longer growth of the 
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disturbances in the computation that results into an earlier development of the roll-over of positive and 
negative disturbance components. 



Experimental data at xj / c=0.60 (Bertolotti, 1996 ) 


Xi jc— 0.50 





XjW 


Xj / C=0.70 [mm] 2 



5 10 15 


Xj [mm] 

a) Steady Case 



30 35 


b) Interaction Case at T— 0 


Fig. 4.11. Total uj -disturbance quantity from a PSE computation 


In order to further describe the purely stationary vortex that will be investigated for its instability to 
three-dimensional high-frequency secondary disturbances in Part 2 of the paper, we show the shape functions 
of the individual modes and the total velocity component along the stationary vortex at different spanwise 
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positions in Figures 4.12 and 4.13, respectively. Here, and w, are the velocity components perpendicular 
and tangential to the constant phase lines (see Figure 2.1). 




Fig, 4.12, Shape functions of individual modes at two streamwise locations 


For the two streamwise positions we selected, it can be seen from Figure 4.10 that the stationary 
disturbances are about to saturate at z, /c=0.6 where they reach amplitudes of A, <u ~ 20%. At z, /c=0.7, 
they are fully saturated at an amplitude level of A, tU ~ 30%. As observed for Swept Hiemenz Flow in Section 
3.2., Figure 4.13 reveals the presence of highly inflectional profiles at both streamwise locations. Comparing 
the results presented in Figure 4.12 with the shape functions for Swept Hiemenz Flow in Figure 3.4, we 
note the very similar shape and even magnitude of the individual shape functions in this region of nonlinear 
saturation of the stationary crossflow vortex. 
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Summarizing the previous plots and visualizing the spatial development of the flow in the region of 
nonlinear saturation, we plot the development of the crossflow vortices for both the purely stationary and 
the interaction case in Figures 4.14 and 4.15. Shown are the contours of the total velocity component in the 
chord wise direction in the region between 50% and 70% chord length, where the nonlinear interaction of the 
disturbance components dominates the physics of the flow. Clearly observed are the developing distortion 
of the boundary layer at x t /c= 0.50, and the roll-up of the vortices in the positive spanwise direction in the 
purely stationary case. For the interaction case at time T— 0, we first note that the disturbances do not 
reach as far into the outer flow and the lift-up of low-speed fluid away from the wall is not as strong as in the 
stationary case. Second, we see a strong vortex core away from the wall at a location where the stationary 
vortex rolls over. 
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Fig. 4.15. Evolution of crossflow vortices : Interaction case at T=Q 


5. Conclusion. In the present paper, we studied the linear and nonlinear stability of crossflow domi- 
nated three-dimensional boundary layers using the PSE. We investigated Swept Hiemenz Flow and the DLR 
Transition Experiment. We described the structure of the stationary vortex in the region of nonlinear satu- 
ration in detail. Investigating the DLR Experiment, we computed the mean flow directly from the measured 
pressure distribution. The linear results for wave angle, wavelength and phase velocity obtained by spatial 
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theory agree reasonably well with the experimental results. Better agreement was found for the disturbance 
eigenfunctions. Prom our nonlinear computations, we could confirm that including a small amplitude travel- 
ing mode in the analysis decreases the saturation amplitude of the physical velocity component significantly, 
and our results match the spanwise averaged experimental data well. For both investigated problems, we 
documented the existence of strongly inflectional profiles along the direction of the statio nary vortex in 
the region of nonlinear saturation. It seems that a very strong fundamental mode and a large mean flow 
distortion are the key-players in setting up the flow for a high-frequency secondary instability that will be 
investigated in Part 2 of this paper. 
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